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ABSTRACT 

O 

o 

^ I We present high resolution imaging observations of a sample of previously 

H ' unidentified far-infrared galaxies at 2; < 0.3. The objects were selected by cross- 

correlating the IRAS Faint Source Catalog with the VLA FIRST catalog and 
the HST Guide Star Catalog to allow for adaptive optics observations. We found 



in 



o 



two new ULIGs (with Lpjji > IO^^Lq) and 19 new LIGs (with Lpjji > IO^^Lq 



&)■ 

Twenty of the galaxies in the sample were imaged with either the Lick or Keck 
adaptive optics systems in H or K'. Galaxy morphologies were determined using 
^ ■ the two dimensional fitting program GALFIT and the residuals examined to look 

' for interesting structure. The morphologies reveal that at least 30% are involved 

■ in tidal interactions, with 20% being clear mergers. An additional 50% show signs 

of possible interaction. Line ratios were used to determine powering mechanism; 
of the 17 objects in the sample showing clear emission lines - four are active 



CL|' galactic nuclei and seven are starburst galaxies. The rest exhibit a combination 

O ' of both phenomena. 

^ • Subject headings: galaxies: infrared — galaxies: interactions — galaxies: evolu- 



tion — galaxies: fundamental parameters — instrumentation: adaptive optics 



1. Introduction 

Luminous infrared galaxies (LIGs; Ljr > lO^^L©), are the best candidates for a link 
between the more powerful ultra-luminous infrared galaxies (ULIGs; Lm > lO^^L©) and 
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normal quiescent elliptical galaxies (Genzel et al. 2001). ULIGs are often interpreted as 
powerful mergers of gas rich spiral galaxies. Good evidence exists for a correlation between 
IR luminosity and the fraction of galaxies which are interacting (e.g., Sanders & Mirabel 
1996). At < IO^^Lq most IR galaxies are single, gas-rich galaxies powered by normal 
star formation, while at Lm > IO^^Lq there is a large increase in the fraction of strongly 
interacting or merging galaxies and an increase in the fraction of AGN-powered galaxies. 
Results of numerical simulations by Mihos & Hernquist (1994) describe the evolution of 
global star formation rate (SFR) for merging pairs of disk + bulge + halo galaxies. At the 
first close approach, star formation is shown to increase slightly over normal levels. When 
the galaxies finally collide, gas is driven into the compact center of the remnant galaxy and 
the SFR increases rapidly to a peak as much as 70 times the initial rate. 

A population of ULIGs at high z have been shown to share some of the properties of 
the local population by studies with Infrared Space Observatory (ISO) (e.g., Sanders 2002, 

and references therein) and with the SCUBA camera on the James Clerk Maxwell Telescope 
(JCMT) (e.g., Bargcr, Cowic & Sanders 1999, and references therein). If that is the case, 
then ULIGs probably played an important role in the star formation history of the universe. 

Several high resolution imaging studies of ULIGs and LIGs have been conducted with 

the Hubble Space Telescope (HST). In particular. Borne et al. (1999) observed 120, z < 0.2 
ULIGs in snapshot mode using the F814W /-band filter and found that virtually all of the 
objects in their sample are interacting or merging, and that as many as 20% contain multiple 
nuclei or arc dense groupings of interacting (soon-to-mcrgc) galaxies. Farrah ct al. (2001) 
observed a sample of 23 ULIGs with the HST WFPC2 camera in V band, and found 87% 
to be interacting. Borne et al. argue, as have others, that there may be an evolutionary 
progression from compact galaxy groups to galaxy pairs to ULIGs to elliptical galaxies. 
These HST data also reveal unresolved nuclei, probably AGN, in 15% of the objects, in 
good agreement with optical and FIR spectroscopic classifications. 

More recently ULIGs and LIGs have been studied with ground based telescopes. Veilleux, 
Kim & Sanders (2002) observed a sample of 118 ULIGs with the University of Hawaii 2.2m 
telescope in R and K'. Optical spectroscopy for this sample was published in Veilleux et al. 
(1999). They find virtually 100% of the sample to be interacting, 39% to be in the early 
stages of merging and 56% to harbor a single disturbed nucleus in the late stages of a merger. 
5% were found to be multiple mergers. They find 35% of their surface brightness profiles to 
be fit well by a pure de Vaucouleurs R^^^ profile and another 38% to be fit equally well by 
either an exponential or de Vaucouleurs profile. Mean half-light radius for their ULIGs was 
found to be ~3.5 kpc in K'. 

Considerable observational effort has been directed at determining whether ULIGs and 
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LIGs are powered mainly by starbursts or active galactic nuclei (AGN), what the evidence 
is for morphological evolution, and whether these correlate with far-infrared (FIR) luminos- 
ity or spectroscopic classification (starburst, AGN, LINER). Recent spectroscopic surveys 
have shown that most of the FIR galaxies seem to be powered by starbursts, but that the 
fraction of AGN-powered galaxies increases with FIR luminosity (e.g., Veilleux et al. 1999). 
A morphological merger sequence which correlates with these spectroscopic classifications is 
not clear, most likely because of large differences in the time scales for the various events. 
Arribas et al. (2004) observed 30 LIGs with the Nordic Optical Telescope in the visible bands 
B, V and /. They find that the LIG population is dominated by starbursts while a higher 
proportion of ULIGs are dominated by AGN activity and could actually evolve into QSOs. 

In an effort to construct larger samples of LIGs and ULIGs having high resolution 
imaging in the near IR, we have identified a new set of candidate objects following the 
method of Stanford et al. (2000) as described below. In addition, we have cross-correlated 
this sample with a bright star catalog so as to allow for adaptive optics (AO) observations 
of the sample. The main aim of our study was twofold: 1) to identify LIGs and ULIGs at 
higher {z > 0.1) redshifts than current FIR-selected samples, and 2) to perform a detailed 
high resolution morphological study that would allow us to identify morphological sequences 
and characterize galaxy interactions in these objects. 



2. Sample Selection 

We have constructed a sample by cross-correlating the FIRST catalog (S'i.4ghz > 1 
mJy, 5(7 with 5" resolution; Becker et al. 1995) with the IRAS Faint Source Catalog {FSC, 
'S'eo/im > 0.2 Jy, 5a; Moshir et al. 1992). The sky coverage for the FIRST catalog is roughly 
from RAs 8 to 17 hrs. 

We extracted all sources which were optically faint, as detailed in Stanford et al. (2000). 
For a flux-flux plot and a plot of radio power at 1.4 GHz versus FIR luminosity illustrating 
the entire cross-correlated FIRST-FSC sample see flgures 1 and 4 in Stanford et al. (2000). 
A major advantage of choosing these FIRST-FSC (FF ) matches is that it also provides a 
good reason to believe that the FIR fiux comes from the optical object at the radio source 
position within the large IRAS error eUipse. 

The FF sample was further cross-correlated with the HST Guide Star Catalog (GSC) 
to define a sub-sample of LIG/ULIG candidates within one arcminute of stars of magnitudes 
brighter than 13 in R. The nearby stars can be used as guide stars for observations with the 
Lick and Keck AO systems. This yielded a sample of ~ 100 targets. 
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Since the aim of this study was to identify new LIGs and ULIGs, we only considered 
those objects for which no pubhshed redshift was available at the time (although the redshift 
of roughly half of the objects in the sample has been published since by the Sloan Digital 
Sky Survey and other authors). The cross-correlation of the three catalogs yielded a sample 
of ~ 50 targets with RAs between 8 and 17 hrs that were previously unidentified. Due to 
observing constraints and weather conditions, we were only able to obtain redshifts for 28 
objects. Of these, two are found to be ULIGs, 19 are LIGs, and the remaining seven we 
simply designate as IR galaxies (IRGs) with FIR luminosities lO^'^ < LpiR < IO^^Lq. 

Finally, wc obtained Lick and Keck AO images of 20 of these objects as described below. 
The sample is given in Table 1. The 20 objects imaged at Lick and Keck observatories appear 
first in the table, in order of increasing RA. The additional 8 galaxies, for which we have 
only spectra, are listed below the horizontal line. 

Column 1 is the target galaxy with its FIRST-FSC catalog name. Columns 2 and 3 
hst, respectively, the J2000.0 RA and DEC of the target. Column 4 is the redshift obtained 
from the object spectrum. Column 5 is the 60//m flux. Column 6 is the 100/im flux; 
values in parentheses indicate upper limits. The fluxes are the template amplitudes from the 
"1002" median scans obtained with the SCANPI utility at IPAC. Column 7 is the luminosity 
distance, calculated using h = 0.71, = 0.73 and Qm = 0.27 (which we assume throughout 
the paper). Column 8 is the integrated flux at 1.4 GHz. Column 9 gives the FIR luminosity, 
calculated as in Stanford et al. (2000). ULIGs and LIGs are normally defined according to 
LjFt (8-1000/im) as a whole. We have based our ULIG and LIG definitions on Lpj^i because 
in the majority of cases we only have IRAS detections at 60/im and lOO/xm for the objects 
in our sample. Most of the objects in the sample were not firmly detected by IRAS at 12//m 
or 25//m. As a consequence it is possible that some of the objects we have classified as LIGs 
may actually be ULIGs. 

The definition for LpiR used is from Sanders & Mirabel (1996) and takes the form 



L(40 - 500/xm) = 4:nDlCFFiR[LQ] 



(1) 



where Dl is the luminosity distance in Mpc, 



Ffir = 1.26 X 10-^^(2.58 x /eo + fioo)[Wm-^] 



(2) 



and C = 1.6. 
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Table 1. The Sample^ 



Object 


Oi 


5 


z 


J GOfJ-TTl 


f b 

J lKj\JfJ,7n 


DL 






Name 


fJ2000~) 


fJ2000) 




(Jv) 


V-'J / 


fMnc) 








(2) 


(3) 
V-') 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


FF0819+2707 


08 19 16.8 


27 07 34 


0.2613 


0.74 


1.07 


1314 


5.97 


12.50 


FF0825+5216 


08 25 34.5 


52 16 42 


0.1726 


0.73 


(0.65) 


823.6 


2.4 


(12.03) 


FF0835+3142 


08 35 51.6 


31 42 00 


0.0483 


0.36 


(0.30) 


211.7 


1.02 


(10.53) 


FF0839+3626 


08 39 50.5 


36 26 57 


0.0961 


0.48 


1.14 


435.9 


3.12 


11.44 


FF0841+3557 


08 41 30.9 


35 57 45 


0.0502 


0.57 


1.68 


220.3 


2.72 


10.98 


FF0934+4706 


09 34 04.0 


47 06 02 


0.1207 


0.33 


(0.58) 


556.8 


1.46 


(11.43) 


FFlllO+3130 


11 10 02.2 


31 30 02 


0.1171 


0.29 


(0.20) 


538.9 


1.06 


(11.23) 


FF1113+5524 


11 13 38.6 


55 24 41 


0.0382 


0.71 


1.50 


166.1 


3.33 


10.75 


FF1122+4315 


11 22 03.6 


43 15 56 


0.1463 


0.46 


0.93 


686.4 


3.47 


11.79 


FFl 138+4405 


11 38 35.5 


44 05 28 


0.0359 


0.68 


0.86 


155.9 


3.35 


10.59 


FF1316+2511 


13 16 42.1 


25 11 56 


0.1459 


0.34 


(0.59) 


684.4 


1.87 


(11.63) 


FF1318+3250 


13 18 24.3 


32 50 41 


0.0367 


0.51 


0.84 


159.4 


2.68 


10.53 


FF1412+4355 


14 12 29.6 


43 55 55 


0.1332 


0.59 


0.85 


619.6 


1.68 


11.75 


FF1429+3146 


14 29 56.5 


31 46 02 


0.1761 


0.18 


(0.25) 


842.1 


1.09 


(11.50) 


FF1517+2800 


15 17 52.8 


28 00 50 


0.1016 


0.39 


(0.60) 


462.6 


2.24 


(11.32) 


FF1519+3520 


15 19 58.4 


35 20 37 


0.1098 


0.24 


(0.49) 


502.8 


2.1 


(11.23) 


FF1708+4630 


17 08 54.0 


46 30 46 


0.2630 


0.29 


0.90 


1323 


2.99 


12.26 


FF1709+5220 


17 09 00.8 


52 20 03 


0.1689 


0.20 


(0.52) 


804.0 


1.18 


(11.61) 


FF1712+3205 


17 12 07.9 


32 05 33 


0.0372 


0.27 


(0.66) 


161.7 


1.19 


(10.34) 


FF1725+4559 


17 25 00.3 


45 59 43 


0.0625 


1.10 


1.03 


276.8 


4.43 


11.26 


FF0834+4831 


08 34 46.8 


48 31 39 


0.1735 


0.30 


(0.80) 


828.3 


2.13 


(11.84) 


FF1439+3232'= 


14 39 16.9 


32 32 39 


0.2502 






1250 


1.21 




FF1601+4514 


16 01 56.6 


45 14 03 


0.0969 


0.58 


1.12 


439.8 


3.35 


11.50 


FF1621+2214 


16 21 08.1 


22 14 08 


0.0843 


0.29 


0.60 


379.2 


1.37 


11.08 


FF1651+3001 


16 51 22.6 


30 01 04 


0.0592 


0.55 


0.96 


275.4 


4.32 


11.02 


FF1656+2644 


16 56 46.5 


26 44 57 


0.1193 


0.44 


(0.69) 


574.8 


1.16 


(11.55) 


FF1721+2951 


17 21 43.7 


29 50 59 


0.1052 


0.23 


(0.48) 


480.2 


1.07 


(11.19) 


FF1723+3845 


17 23 29.7 


38 45 12 


0.0377 


0.22 


(0.28) 


163.9 


1.04 


(10.15) 
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^Horizontal line divides imaged and non-imaged samples. 
''Parentheses indicate upper limits. 

'^IRAS data for this object suffers contamination from a nearby M6 star. 
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3. Observations and Data Reduction 

Spectroscopic observations of each galaxy in the sample were obtained using the Kast 
Double Spectrograph at the cassegrain focus of the Shane 3-meter telescope at Lick Ob- 
servatory. We used the 600/4310 grism for the blue side and the 600/7500 grating for the 
red side to obtain a useful wavelength coverage spanning from the atmospheric cutoff at ~ 
3400 A to 8100 A. We used different sht widths to match the seeing conditions, typically 
between 1" and 2'/5, yielding a resolution between 2.6 and 6.5 A pixel"^ for the blue side and 
between 3.3 and 8.3 A pixel"^ on the red side, so that the typical resolution for the spectra 
was roughly 300 km s~^. The total integration time for each galaxy was 900 s. 

The spectra were reduced with IRAF, using standard reduction procedures. After cor- 
recting for bias and fiat fielding, we subtracted the sky and wavelength-calibrated the two- 
dimensional spectrum using OH skylines (for the red side) and arc lamps (for the blue side). 
We then flux-calibrated the spectra using spectrophotometric standards from Massey et al. 
(1988) and extracted the spectra using the IRAF apextract routines. In the majority of 
cases, we measured redshifts from stellar absorption lines, so that the redshifts we quote 
correspond to the stellar component of the galaxies as opposed to the gas. 

The galaxies FF 1122+4315 and FF 1429+3146 were imaged in K' using the Keck II AO 
system (Wizinowich et al. 2000a,b; Johansson et al. 2000) with the NIRC-2 camera (PI: K. 
Matthews & T. Soifer). Both galaxies were observed using the NIRC-2 Wide-Field camera, 
which yields a plate scale of 0'.'04 pixel"^. In addition, FF 1429+3146 was also observed with 
the Narrow-Field camera, which yields a plate scale of O'.'Ol pixel"^. The remaining galaxies 
in the sample were observed using the natural guide star LLNL AO system on the 3 meter 
Shane telescope at Lick Observatory; for details about the LLNL AO system refer to Bauman 
et al. (1999) and Gavel et al. (2000). The AO system feeds the AO-optimized infrared camera 
IRCAL (Lloyd et al. 2000), yielding a plate scale of 0'.'076 pixel"^. Observations at Lick were 
done in H rather than Kg since the warm optical elements in the AO system result in a high 
thermal background in the latter. The AO FOV is 20" for the Lick system, 10" for the Keck 
Narrow-Field camera and 40" for the Keck Wide-Field camera . 

Observations of point spread function (PSF) stars were obtained either immediately 
before or after the observations for each galaxy. We attempt to account for anisoplanatism 
by matching the distance and position angle from the guide star (GS) to the PSF. However, 
atmospheric conditions vary somewhat on a shorter scale than our total integration times, so 
that each PSF is close, but does not perfectly match the conditions for each image. For this 
reason, we are unable to provide precise Strehl ratios for each image, but we estimate that 
the typical ratios for all our images were between 0.1 and 0.2. We provide full width at half 
maximum (FWHM) of the PSF for each image as an indication of the system performance 
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for every field. Only the observations at Keck and those at Lick on 25 Jan 2003 were 
done under photometric conditions; we estimate an extinction in H between 0.1 and 0.4 
magnitudes for the rest of the observations. The images were reduced with IRAF, using 
standard IR reduction procedures. A complete journal of observations is given in Table 2 
which includes total exposure times and GS information. 



Table 2. Journal of Observations 



Object 


Scale (kpc/'O 


GS V 


Separation (") 


PA (deg) 


Exp. Time (s) 


FWHM (") 


Date Obs. 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


FF0819+2707 


4.028 


11.5 


17.9 


64.8 


12x300 


0.49 


25 Jan 2003 


FF0825+5216 


2.904 


11.9 


30.2 


61.2 


14x300 


0.49 


25 Jan 2003 


FF0835+3142 


0.981 


12.5 


37.9 


286.1 


10x300 


0.49 


25 Jan 2003 


FF0839+3626 


1.809 


11.7 


31.3 


250.6 


6x300 


0.15 


03 Apr 2002 


FF0841+3557 


0.968 


12.5 


31.6 


151.6 


3x300 


0.49 


25 Jan 2003 


FF0934+4706 


2.149 


12.8 


34.8 


5.4 


12x300 


0.31 


03 Apr 2002 


FFlllO+3130 


2.134 


12.6 


43.6 


127.2 


6x600 


0.80 


19 Mar 2003 


FF1113+5524 


0.762 


7.5 


37.0 


295.9 


7x300 


0.49 


25 Jan 2003 


FF1122+4315'^ 


2.563 


12.7 


27.3 


267.9 


25x60 


0.068 


25 May 2002 


FF1138+4405 


0.732 


9.7 


29.6 


258.9 


10x300 


0.43 


03 Apr 2002 


FF1316+2511 


2.527 


11.2 


40.6 


79.9 


8x300 


0.60 


03 Apr 2002 


FF1318+3250 


0.719 


10.5 


34.0 


115.5 


6x300 


0.49 


03 Apr 2002 


FF1412+4355 


2.371 


9.8 


27.4 


350.7 


10x300 


0.39 


03 Apr 2002 


FF1429+3146^ 


3.012 


11.9 


30.6 


8.2 


10x120 


0.092 


25 May 2002 












5x120 


0.068 


25 May 2002 


FF1517+2800 


1.848 


11.9 


30.0 


13.6 


6x300 


0.42 


03 Apr 2002 


FF1519+3520 


1.979 


11.1 


36.2 


188.3 


6x300 


0.29 


19 Mar 2003 


FF1708+4630 


4.024 


9.1 


27.2 


161.2 


7x300 


0.44 


03 Apr 2002 


FF1709+5220 


2.853 


11.8 


27.2 


237.0 


11x300 


0.14 


04 Sep 2004 


FF1712+3205 


0.729 


11.1 


43.0 


254.1 


5x300 


0.17 


14 Aug 2003 


FF1725+4559 


1.189 


11.9 


27.4 


187.5 


16x300 


0.63 


05 Sep 2003 



^Galaxies observed in K' with the Keck II telescope 



- 11 - 



4. Analysis 

4.1. Fitting Technique 

Fitting a mathematical model to an image of a galaxy is one consistent way of determin- 
ing its morphology Several different mathematical models have been used over the years to 
fit the most common galactic shapes. These include the well known de Vaucouleurs profile 
that models elliptical galaxies and the exponential profile for galactic disks. More recently 
the Sersic (1968) model has become a highly valuable tool for modeling various components 
of galaxies including bulges, disks, and bars. Using the Sersic model is beneficial because 
the model is able to adapt to the de Vaucouleius (elliptical) profile at Sersic index N=4, the 
exponential disk profile at N=l, and a Gaussian shape at N=0.5. 

Many authors have used a one dimensional ellipse fitting routine to plot a light profile 
such as IRAF ellipse. For example, Veilleux, Kim & Sanders (2002) use the standard Fourier 
expansion of Binney & Merrifield (1998) to fit isophote ellipses to their galaxies. After 
generating these ellipses the programs then make a plot of isophote intensity versus radius 
and derive the surface brightness profiles by making a best fit to those points. Surveys 
based on these routines usually classify the object as either elliptical or disk shaped based 
on whether a de Vaucouleurs or exponential model fits best. While this is not unreasonable, 
some of these classifications now need to be revised because a one dimensional fit can be 
subject to errors due to isophote twists and the large variety of galaxy morphologies. Some 
of these surveys use a profile slice along the major or minor axis and some use both. The 
profile can be different depending on whether the major or minor axis is used (see Peng et 
al. 2002, and references therein). 

The galaxies in this sample were fit in two dimensions using a program called GALFIT 
(Peng et al. 2002). The program bypasses ellipse fitting, and fits a model to the light profile 
directly. GALFIT uses fitting to minimize the error in two dimensions and can handle the 
fitting of multiple models simultaneously. The best fit is then subtracted from the data, and 
the residuals can be analyzed. In this survey, one or two models were used as needed to make 
a reasonable subtraction. More models can be used, but improvements to the subtraction 
are not always desirable. Generally, as resolution improves, many Sersic profiles can be fit 
to a galaxy at the same time but these are not necessarily physically meaningful. Depending 
on the number of components, a Sersic model can have anywhere from eight to several dozen 
free parameters. 

The Sersic parameter N affects the degree of cuspiness of the galaxy. As the Sersic index 
decreases, the galaxy becomes more "cuspy" in the center. This means the core intensity 
flattens quickly as r increases to the half-light radius and the intensity falls off steeply beyond 
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the half-light radius (Peng et al. 2002). Unperturbed spiral galaxies are a homogeneous 
group. Their Sersic index usually does not deviate from N=l. Similarly, ellipticals do not 
deviate much from N=4, though they are a less homogeneous group. Mergers of two disk 
galaxies tend to produce a bulge with knots of compact unresolved star formation and some 
mass concentrated in the center. In general, these types of mergers can have Sersic indices 
that are greater than 4 or other ambiguous profiles. 

4.2. Fitting Procedure 

The fitting process consists of two stages. The first stage involves an initial guess. A 
visual evaluation of the galaxy is made to decide what types of objects arc present (disks, 
bulges, bars, etc). The general parameters of these objects, such as effective radius, centroid, 
ellipticity and position angle can all be approximated by examination of the image with an 
image analysis tool. 

The fitting program uses the initial guess parameters to create a model, and convolves it 
with a PSF provided by the user in order to match the seeing conditions and the resolution. A 
"best fit" model is obtained by performing a least squares fit of the PSF-convolved model to 
the data. The second stage consists of refinements to the model to achieve the most accurate 
representation of the galaxy. The accuracy of the model can be assessed by examining the 
residual image for artifacts of under- or over-subtraction. The residual image is the image 
formed when the model is subtracted from the original image. If the galaxy is simple and 
unperturbed then the residuals should be only the sky background or spiral arms for spiral 
galaxies. This is especially true for distant galaxies where the resolution that we achieve 
prevents us from resolving extreme detail such as globular clusters. 

If the galaxy is perturbed then other types of residuals will be apparent. These may 
consist of dust lanes, tidal tails, or multiple nuclei. Large symmetric areas with negative 
values arc characteristic of over-subtraction. These types of errors are generally easy to spot 
and another fit should be considered. Sometimes the model may be obviously wrong. Many 
adjustments to the input parameters can be made to refine the model. One possibility is 
isolating variables by fixing quantities that are known to be accurate and letting the program 
fit a smaller subset of the parameters. 

The accuracy of the parameters determined by the fitting procedure will naturally de- 
pend on how well the PSF matches the conditions under which a given image was taken. 
To go beyond a simple determination of disk versus elliptical and look at more complicated 
features, we need to evaluate the uncertainty in the model due to the PSF. We ran several 
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tests where we used a "mismatched" PSF to fit various galaxies. The mismatched PSFs were 
obtained with the same instrument and set up, but under different observing conditions and 
with FWHM that differed by as much as 50% of the value of the actual corresponding PSF. 
We found that, overall, the resulting model parameters changed by only a small percentage 
of their value. In particular, Sersic indices of N > 1 changed by less than 5% and the position 
angle by only a few degrees. Although smaller Sersic indices (N < 0.5) varied by a larger 
percentage (as much as 40%), their values remained close to the values for disk profiles. The 
effects are similar to running the program with no PSF. It should be noted that in some cases 
the PSF can be so "bad" , for instance when position angle differs greatly from the image, 
that GALFIT will not be able to converge. In such cases it is better to run the program with 
no input PSF. Wc therefore feel confident that, while the precise value of the Sersic index 
may be uncertain due to uncertainties in the PSF, the overall determination of the galaxy 
morphology and global features is robust. A more serious concern when a PSF is mismatched 
is that the structure in the residuals will be affected. If the PSF does not properly represent 
the image quality of the galaxy, fine structure in the residuals will probably be lost. 

Another factor that will affect the Sersic parameter is surface brightness fading. For 
galaxies at higher the disk of a galaxy (where there is less signal) will fade faster than the 
central regions causing the galaxy to appear significantly smaller. This is because flux scales 
roughly as where L is the luminosity measured at the source and z is redshift (Misner, 

Thorne & Wheeler 1973). In the case where there is a disk + bulge, losing the edges of the 
disk will increase the cuspiness and raise the Sersic index. 

Finally a limited fleld of view may effectively raise the Sersic index. When combined 
with surface brightness fading, a situation is created where a galaxy appears to be smaller 
and the drop off from the center is steeper. In such a case only the central bulge of a galaxy 
is modeled. For nearby galaxies where the field of view is less than the scale length r^, the 
model for the galaxy may not be correct (Peng et al. 2002). Even for our nearest galaxy, 
FF 0841+3557, at ^ = 0.036 the FOV for the Lick AO system corresponds to 14 kpc. For the 

exponential profile Vg — (1.678r 1 ), and if we take FF 0841+3557 as an example, we calculate 

2 

Tg to be 7.7 kpc. 

Figure 1 is an example of how a residual image is obtained for two galaxies in the sample. 
The top three panels, showing FF 1110+3130, are an example of a numerical model that is 
an accurate representation of the morphology of the galaxy. It produces a clean residual 
which shows sky background. Panel a is an image of the galaxy before subtraction. Panel 
6 is a model of the galaxy and panel c is the residual image produced from the subtraction 
of panel h from a. The bottom three panels, showing FF 1519+3520, are an example of a 
reasonable fit with interesting residuals, namely a compact companion ~ 0'.'5 east of the 



-14- 



galaxy nucleus. In this case we have only modeled the larger galaxy to the west. 




Fig. 1. — For both rows, panel a is an image of the galaxy before subtraction. Panel b 
is a numerical model of the galaxy and panel c is the residual image produced from the 
subtraction of panel b from a. In this and the following figures, north is up, east is to the 
left, and the scale bar represents approximately 1 kpc. 



4.3. Results 

Figure 2 shows the H or K'-hand image of each of the 20 galaxies imaged in the sample. 
In terms of FIR luminosity, the 20 imaged objects consist of two ULIGs, 13 LIGs and five 
IRGs. The additional eight galaxies for which we only have spectra yielded six LIGs and 
one IRG. 

For each of the galaxies modeled in the sample, the parameters for their models are 
listed in Table 3. Column 1 lists the target galaxy which may have one or two component 
objects that are modeled within it. Column 2 lists the apparent H or K' magnitude of the 
object being modeled as determined by the best fit. Column 3 is the Sersic parameter, N. 
When the model used was an exponential profile the value of the Sersic index is by default 
set to 1. The Sersic model or "hght profile" is given by: 

-K[(r/ri)l/iV_i] 

E(r) = EeC 5 (3) 

where E(r) is the surface brightness at a given radius r, ri is the effective radius of the 
galaxy. Eg is the surface brightness at r i , N is the Sersic index and k, is coupled to N so that 
half of the total fiux is always within ri. This model is described in detail in Peng et al. 
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(2002). The paper should be consulted for the mechanics of how the model works, because 
more equations are involved than what is shown here. The exponential profile is given by: 

S(r) = Soe-(''/'^^) (4) 

where Tg is the scale length (1.678ri). Column 4 lists the half-light radius (ri) of the object 
in kpc. Column 5 is the ellipticity {e) of the isophotes (defined as the 1 — b/a where b/a is 
the axis ratio for an ellipse). Column 6 is the position angle in degrees. Column 7 is the 
boxiness-diskiness parameter. A negative value corresponds to a disky galaxy and a positive 
value corresponds to a boxy galaxy (see Section 4.4 below). Column 8 is the model type used 
for the object where "Ser" is short for Sersic and "Exp" is short for exponential. Column 
9 indicates whether the galaxy appears to contain a single nucleus, a double nucleus, or 
multiple nuclei. 

Column 10 indicates whether the galaxy is interacting or not. We classify a galaxy as 
interacting ("Y") if the galaxy appears to be in the early stages of a merger (i.e., still showing 
distinct components) and/or shows obvious tidal tails or other tidal debris. We designate a 
galaxy as possibly interacting ("?") if the galaxy was difficult to fit, implying that it may 
have an irregular morphology. The rest we designate as non-interacting ("N"). The models 
used for each individual object are described in detail in Section 5 below. 



4.4. Powering Mechanism 

ULIGs and LIGs can be classified in two categories according to the primary source of 
their luminosities: those galaxies that achieve their high luminosity from starburst activity 
(henceforth referred to as a starburst galaxies) and those that are mainly powered by AGN. 
It is plausible that most LIGs contain some combination of excitation mechanisms including 
AGN, starbursts, shocks, mergers, and bars. The question then becomes which mechanism 
is dominant in each galaxy and whether there are any trends evident in the sample. 

Figure 3 shows the emission line flux ratio [0III]/H/3 versus the ratio [Nil] /Ha plotted 
for the 17 objects in the sample which had firm detections of all four emission lines and 
a redshift z < 0.2 (so that both Ha and H/9 are in our observed spectral range). This 
plot is similar to the plot in Fig. 1 of Kauffmann et al. (2003) which itself is derived from 
the "BPT diagram". Baldwin, Phillips & Terlevich (1981) demonstrated that it is possible 
to distinguish type-2 AGNs from normal star forming galaxies by plotting their emission 
line ratios and this idea was expanded upon by others. The curved line represents the 
demarcation between starburst- and AGN-powered galaxies as determined by Kewley et al. 
(2001). 
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According to the diagram, the galaxies FF 0834+4831, FF 0835+3142, FF 1122+4315, 

and FF 1519+3520 arc identified as AGN, while FF 0934+4706, FF 1110+3130, FF 1656+2644, 
FF 1709+5220, FF 1712+3205, and FF 1318+3250 are identified as starburst galaxies. Ad- 
ditionally, FF 1138+4405 lies near the demarcation line, and is likely a star forming galaxy 
as well. In both of these groups, there is a range of FIR luminosities so that we do not 
see a clear correlation between Lpm and powering mechanism for the LIGs and IRGs in 
the sample, although it is important to note that some of the objects are classified as LIGs 
based on upper limits to LpiR and this adds scatter to any possible trends. Neither of the 
two ULIGs are plotted in Fig. 3 because [N II] and Ha were redshifted out of our observed 
spectra. However, they both have [0III]/H/9 ratios characteristic of AGN. Studies with 
larger samples like Veilleux et al. (1999) and Arribas et al. (2004) indicate the fraction of 
AGN dominated LIGs increases with Lpjfi. FF 1651+3001 is the only galaxy that falls in 
the LINER region. The rest of the objects fall in a region intermediate between starburst 
and AGN, possibly representing a population containing some combination of both AGN 
and starbursts, with neither being clearly dominant. 



4.5. Morphologies 

One of the two ULIGs, FF 1708+4630, is a merger at an early stage, while the other one, 
FF 0819+2707, has possible signs of interaction. Three LIGs, FF 1412+4355, FF 1429+3146 
and FF 1519+3520, are also mergers at an early stage, where the two nuclei are still distinct. 
The LIG mergers are likely mergers of two disks since both components can be fit well 
by near exponential profiles and are in early merging stages, because both nuclei are still 
distinct. Between 30 and 60% of LIGs are reported to be mergers in the literature (Sanders 
& Mirabel 1996). Similarly, 23% (3/13) of the LIGs in our sample are found to be mergers. 
In agreement with the rarity of multiple mergers cited by Veilleux, Kim & Sanders (2002), 
we find only one object FF1429+3146 (a LIG) that appears to be a multiple merger in the 
residual image. None of the 5 IRGs have multiple nuclei or obvious signs of interaction. 

There are only 2 objects that seem to have large bulge components in the sample. One of 

these objects is the merger ULIG FF 1708+4630 with an unusually high Sersic parameter at 
N=6. Another bulge dominated object is the IRG FF 0835+3142. Otherwise, the majority 
of the objects seem to be disk-dominated rather than bulge dominated, though the Sersic 
parameters range from 0.7 to 2.3. 

The statistical breakdown for profiles produced in Veilleux, Kim & Sanders (2002) for 
single nucleus objects are that 2% are best fit by a pure exponential disk, 35% are best fit by 
a pure elliptical, and 38% are fit equally well by both. There are 10 single-nucleus LIGs in our 
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Table 3. Model Parameters & Morphologies 



Object 


MAG 


N 


n (kpc) 

2 




PA (deg) 


Boxy/Disky 


Fit Type 


Nucleus 


Interact? 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


f U819+27U7 


25.1 


2.7 


5.5 


A CO 

U.52 


''70 A 

76.0 


~ 


ber 


sinj 


?le 


( 


i" U8z5+5z1d 


25.U 


2.3 


A AA 

0.9U 


C\ A A 

U.44 


—64.2 


A C\>-f 

0.27 


ber 


sinj 


^le 


( 




23.3 


4.0 


0.8 


A OA 

U.80 


1 A 

—61.0 


A OQ 

0.68 


her 


sinj 


;le 


o 
f 


FF0839+3626 


20.0 


1.8 


3.1 


0.52 


^ A A 

—34.4 


—0.31 


Ser 


sinj 


de 


AT 


FF0841+3557 


23.6 


1.0 


4.6 


0.26 


A A 

9.0 


—0.83 


Exp 


sin^ 


de 


TVT 

JN 


1 1 Uy34+47(Jo 


19.0 


7.2 


O A 

3.0 


0.74 


OA O 

82.8 


—0.54 


ber 


sinj 


de 


( 


FFlllO+3130 


22.4 


0.7 


2.9 


0.78 


23.3 


—0.15 


Ser 


sinj 


de 


o 

f 


rb 1113+5524 


21.8 


5.0 


A O 

2.3 


A /IT 

0.47 


17.3 


Air" 

0.15 


Ser 


sinj 


de 


TVT 

JN 




30.5 


1.0 


8.7 


0.64 


O A 

8.0 


—0.81 


Exp 








TT^TT^ 1100 1 /loir" 

t b 1122+4315 


8.8 


3.1 


O A 

6.9 


A CT O 

0.53 


13.1 


A 1 

—0.16 


ber 


sinj 


de 


Y 




on c 

20.5 


1.16 


1.7 


r\ o A 

0.84 


34.6 


A AO 

0.02 


he r 


sin^ 


de 


Y 


FF1316+2511 


20.6 


10.0 


2.8 


0.84 


—33.3 


0.36 


Ser 


sinj 


de 




FF1318+3250^ 














Bar 


sinj 


de 


N 


FF1412+4355 


23.0 


1.7 


0.67 


0.90 


-43.0 


0.47 


Ser 


double 


Y 




21.6 


1.0 


1.9 


0.09 


85.9 


~ 


Exp 








FF1429+3146 


15.2 


2.3 


5.9 


0.52 


-1.4 


-0.51 


Ser 


mult 


Y 




19.1 


1.1 


1.2 


0.68 


24.1 


-0.33 


Ser 








FF1517+2800 


20.7 


2.2 


0.90 


0.61 


-22.6 


-0.23 


Ser 


sinj 


de 


Y 




22.0 


0.3 


2.8 


0.49 


44.7 


0.08 


Ser 








FF1519+3520 


20.3 


1.4 


2.8 


0.68 


-35.4 


-0.20 


Ser 


double 


Y 




23.7 


0.9 


0.72 


0.85 


86.7 


-0.33 


Ser 








FF1708+4630 


21.3 


6.0 


4.0 


0.26 


62.8 


0.31 


Ser 


double 


Y 


FF1709+5220 


23.1 


0.7 


2.3 


0.34 


35.7 


-0.25 


Ser 


sinj 


de 


? 


FF1712+3205 


20.4 


2.5 


1.4 


0.25 


-79.9 


-0.95 


Ser 


sinj 


de 


? 


FF1725+4559 


20.8 


1.1 


1.4 


0.89 


30.8 


0.02 


Ser 


sinj 


de 


? 



^The model parameters for this object are not meaningful since only the central bar was imaged. 



-18- 



sample. Four (40%) are fit fit well by near exponential shapes, and four (40%) fall somewhere 
in between exponential and de Vaucouleurs. The "in between" Sersic indices indicate that 
they can be fit by either a pure exponential or a de Vaucouleurs profile but neither is 
ideal. The remaining two objects (20%) have indices that are unusual high, indicating that 
neither an exponential nor a de Vaucouleurs profile would be a good fit. We do not find any 
preference toward elliptical profiles for the LIGs, and instead find more disks, but this may 
be due to small sample statistics. We obtain a mean r i of 2.9 kpc for the single-nucleus LIGs 
in H (except for FFl 122+4315 which was imaged in K'); this is somewhat smaller than the 
mean half-light radius of 3.5 kpc for ULIGs found by Veilleux, Kim & Sanders (2002). Of 
the IRGs only two are fit well by models. FF1113-I-5524 has a classic disk-|-bulge profile and 
FF 1712-1-3205 has a nearly a lenticular shape. 

Boxiness refers to the shape of the isophotes. Objects with a positive value ("boxy" 
objects) have slightly more square isophotes while those that have a negative value ( "disky" 
objects) are rounded. Boxiness tends to be a sign that a galaxy has undergone a recent tidal 
interaction. Most elliptical galaxies are disky; elliptical galaxies that are boxy tend to have 
higher mass to light ratios (Kormendy & Djorgovski 1989). Boxiness does not seem to be 
correlated with infrared luminosity in this sample. 

Galaxies that are interacting make up 30% of the sample, including one ULIG and five 
LIGs. If we include all objects that are possibly interacting, they would make up 80% of the 
sample. Two of the galaxies that do not show signs of interaction arc IRGs and the other 
two arc LIGs. So, in agreement with previous studies, we observe a trend of toward a higher 
interaction rate at higher luminosities. 

FF 0834+4831, FF 0835+3142, FF 1122+4315, and FF 1519+3520 arc identified as AGN. 
FF 1519+3520 is an early merger of two galaxies, and possibly a dust obscured AGN. 
FF 0834+4831 was not imaged. The other two objects, as well as FF 0819+2707 and 
FF 1429+3146, have point-like cores in their residuals that are likely to be their active nuclei. 
Except for FF 0934+4706 which has a PSF like core, and FF 1318+3250 for which we only 
have an image of the central region, those galaxies below or slightly above the demarcation 
line for starburst galaxies in Fig. 3 can be modeled with near exponential profiles, and none 
of them show overt signs of tidal interaction. Conversely, all of the objects that are currently 
involved in a tidal interaction are found either in the Seyfert region or in an intermediate 
region between the starburst and Seyfert regions. It is possible then that every object in the 
sample that is undergoing a merger has some level of nuclear activity. 
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5. Notes on Individual Objects 
5.1. The ULIGS 

FF 0819+2707 This ULIG is the most powerful object imaged in the sample, with 
a FIR luminosity of log(LFiR/LQ) = 12.5. Unfortunately, we do not have a flux value for 
[N II] or Ha for this object since these lines were redshifted out of the spectra we obtained. 
However, the value of log [0III]/H/3 = 1.0 for this object indicates it would almost certainly 
be in the AGN Seyfert range of the BPT diagram. The spectrum, shown in Fig. 4, shows a 
young stellar component. No signs of an ongoing merger arc evident in the residual, nor is a 
close companion apparent. However, the model is not a good fit since it undcrsubtracts at 
the edges, possibly indicating that the galaxy is perturbed. The redshift for this object was 
pubhshed by Brand et al. (2003) after we carried out our observations. Brand et al. indicate 
that the object lies within a superstructure of radio galaxies. 

FF 1708+4630 For this galaxy, a merger is evident in Fig. 5 after subtraction of a 
boxy elliptical host with Sersic parameter N=6. The merger components are difficult to 
model. One of them is probably the core of the elliptical galaxy that is less than 1 kpc in 
diameter (and correspondingly less luminous). The intruder seems to be a smaller galaxy 
about 3 kpc in diameter. It lies only 1.5 kpc (in projection) from the core of the primary 
galaxy. The object has a starburst spectrum with ongoing star formation as indicated by 
strong O II emission, possibly triggered by the merger. Because of the higher redshift for 
this galaxy, [N 11] and Ha are also redshifted out of our spectral range. However, the value 
of log [0111]/H/5 = 0.62 indicates that this object is likely in the AGN region. The Mgl6 
and Ca II absorption features are present indicating that an older population is also present. 

5.2. The Lies 

FF 0825+5216 The dominant power source of this log(Lpif>/LQ) < 12.04 galaxy is 
uncertain. Its spectrum is that of an old population with very weak [0 11] emission. The 
galaxy is less than 10 kpc in diameter in projection and is best modeled by a Sersic model 
with N=2.29. No signs of a companion are observed. The model is fairly accurate and there 
are almost no residuals that are detectable. The galaxy may have gone through a ULIG 
phase in the past and is now relaxing to an elliptical shape. 

FF 0839+3626 After subtraction, this moderately luminous LIG reveals a barred spiral 
(SBc). The disk is modeled by a Sersic profile with N=1.77. Based on the spectrum, the disk 
is probably composed of a population of older stars. The larger arm is about 2.1 kpc wide 
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and 8.4 kpc long. This arm is actually visible before subtraction. The other arm is short 
and faint. The bar is about 1.4 kpc wide. No companion is visible. The spectrum reveals 
low ionization emission lines with a strong Ha line. Unfortunately, [O 111] fell precisely in a 
small gap between the blue and red detectors so that we cannot obtain an [0111]/H/9 ratio. 
However, the value of [NII]/Hq; = —0.16 indicates that the object may fall in the AGN or 
transition region. 

FF 0841+3557 This object lies on the boundary between LlGs and IRGs. The galaxy 
does not seem to be interacting. The spectrum indicates it has an old population. There are 
practically no residuals in the subtracted image. It is best modeled by an exponential profile. 
From the morphology it appears to be an ordinary disk galaxy; no bulge is detectable, but 
this could be due to undersampling. 

FF 0934+4706 The spectrum of this object shows a blue continuum with low ionization 
emission lines. Based on its line ratios, the object has a starburst. This galaxy has a cuspy 
profile reflecting a high concentration of mass in the center. It is difficult to model accurately. 
It also has a bar-like structure visible in the residual but spiral arms are not apparent. 

FF 1110+3130 This lower luminosity LIG shows a starburst spectrum with a blue 
continuum. The galaxy is fit accurately by a Sersic profile with N=0.67, and there are no 
identifiable residuals. The shape of this model is near Gaussian (the Sersic profile assumes the 
Gaussian shape when N=0.5). The galaxy is probably a perturbed disk, but no companion 
galaxy is visible. 

FF 1122+4315 This luminous LIG appears to be a spiral galaxy with a bulge-like 
component that extends over most of the image. The spectrum indicates it is an AGN 
Seyfert galaxy. An elliptical shape with a Sersic index of N=3.09 was used to model the 
bulge component. The residual image (Fig. 6) shows much structure including a possible 
warped disk and/or tidal debris, although some effects of over-subtraction are evident. A 
small knot to the southwest could be the core of an interacting dwarf galaxy or a bright knot 
of star formation. 

FF 1316+2511 This is a strongly perturbed LIG. The Sersic exponent, while unusu- 
ally high (N=10), seems to produce an accurate model of this cuspy galaxy; practically no 
residuals are seen after subtraction. The galaxy is spherical in shape (EG). It is boxy, but no 
other signs of recent interaction are apparent. The spectrum indicates that an old population 
of stars is present. 

FF 1412+4355 This object is a likely starburst galaxy though it hes slightly above the 
demarcation line in Fig. 3. Figure 7 shows that the galaxy is an early merger of two disk 
galaxies with tidal tails and extended debris. The larger spiral galaxy is well modeled by 



- 21 - 



an exponential profile, and it has a very long tidal tail. The galaxies are ~1.8 kpc apart in 
projection. The smaller companion is modeled well by a boxy Sersic with N=1.7, and shows 
practically no residuals after subtraction, while the larger galaxy reveals a distinct core 0.8 
kpc in diameter and a diffuse cloud of debris. 

FF 1429+3146 The emission line ratios for this object indicate it is somewhere between 
the AGN and starburst regions on the BPT diagram. However given the morphology of 
this galaxy, and a possible faint broad emission line component in its optical spectrum, it 
could be an obscured AGN. The galaxy is clearly undergoing an early merger and contains 
a young population of stars. We obtained two separate images of this galaxy using the 
narrow and wide NIRC-2 cameras. Figure 8 shows the narrow camera image. Each image 
was modeled independently; the model parameters determined for each of these images are 
generally consistent with each other, although the (higher resolution) narrow camera image 
yields a better fit. The parameters listed in Table 3 correspond to the narrow camera 
image. FF 1429+3146 has four distinct features. The two major components are a large 
cuspy Sersic profile and a smaller merging companion to the south also with a Sersic profile. 
There is a dwarf companion or star forming region to the northeast and possibly one to the 
west. All of these components are visible before subtraction. The north galaxy has a Sersic 
parameter of N~2, while the south galaxy has a Sersic parameter of N~l.l corresponding 
to an exponential shape. The models for both major components show diskiness. In the 
residual images, a core about 0.5 kpc across corresponding to the larger galaxy can be seen. 
This is possibly an active nucleus. The galaxies are about 1.5 kpc apart in projection. The 
unusual shape of this galaxy could indicate it is undergoing a major merger. In addition 
there may also be a minor merger with a smaller companion making it a multiple merger. 
There seems to be some dust present to the north and south. 

FF 1517+2800 This galaxy has a bright Sersic profile with N=2.17. The galaxy also 
has a faint Sersic component ~0.5 kpc off center with N=0.31. This unusual shape is similar 
to a Gaussian, but steeper. The subtracted image shows a core ~1 kpc in diameter and some 
other residuals that may be a tidal tail. Some debris to the south of the obvious core may 
be another core that is part of a former companion. The spectrum of this galaxy is heavily 
reddened and shows strong emission lines. Emission line ratios indicate that the object is 
probably a starburst galaxy, although it hes in the region between starbursts and AGNs in 
Fig. 3. 

FF 1519+3520 This system is best modeled by disk components. The larger com- 
ponent has a Sersic index of N=1.44 while a smaller component on the southeast has an 
N=0.91 index and is ~4 magnitudes less luminous. There is a bright, resolved core about 
1.5 kpc in diameter roughly at the position of the smaller galaxy. It is presumably merging 
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with the primary source. There are features that could be tidal tails to the north and south. 
The spectrum indicates that the stellar population is older. Its FIR luminosity puts it at 
the lower end of the LIG scale. 

FF 1709+5220 This moderately luminous LIG is likely to be powered by a starburst. 
The galaxy is best modeled by a Sersic profile of index N=0.66; this index could be low due 
to undersampling. The image is not deep enough to determine if there is a bulge component. 
The galaxy may have a companion of similar radius and magnitude, with a Sersic index of 
N=0.48 that lies 16" away in the image. Since we do not have a redshift for this object, we 
cannot determine whether it is a projected galaxy or a true companion. 

FF 1725+4559 This low luminosity LIG is best modeled by an exponential profile. It 

is slightly boxy but otherwise appears normal. The subtracted image reveals a possible core 
about 1 kpc long and 0.5 kpc wide. It lies at the same position angle as the disk and it could 
be a bulge that is too small to model. In that case this galaxy would seem to be a normal 
disk + bulge galaxy. 

5.3. The IRGs 

FF 0835+3142 This object contains an older population of stars. It can be fit by a 
boxy elliptical host but is difficult to model accurately. The subtracted image reveals a 
possible unresolved core. 

FF 1113+5524 This object seems to be a normal spiral galaxy. Its spectrum shows 
some some star forming activity. The residuals show that it is a prime example of an Sbc 
spiral galaxy. A prominent bar and two arms nearly 1 kpc wide each can be seen. 

FF 1138+4405 This object can be fit with a near exponential profile. The residual 
image reveals a very unusual oblong shaped core. The spectrum of this object is reddened 
showing very strong emission (the strongest of the sample) that classifies this object as a 
starburst, and an underlying old stellar population. 

FF 1318+3250 The image of this galaxy in Fig. 2 shows only the central bar of the 
galaxy. The galaxy has two faint spiral arms visible only in optical images. The bulk of the 
NIR fiux comes from the bar, which is ~15 kpc long in projection. The bulge is ~4 kpc in 
diameter. The position angles of the bar and bulge are offset by 20 degrees. The spectrum 
of this galaxy shows some star formation based on the presence of [O II] emission, but the 
bulk of the stellar population is probably old. No AGN is apparent. 

FF 1712+3205 This galaxy can be modeled by an N=2.5 disky Sersic, which is roughly 
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a lenticular shape. The spectrum shows weak Ha and [O II] emission. 

6. Summciry 

The results from this study indicate that the technique of cross-correlating 1.4 GHz and 
far-infrared fluxes selects many perturbed galaxies, with frequent early mergers and merger 
remnants. It selects a high number of starburst galaxies at redshifts 0.1 < 2; < 0.3. A few 
type-2 AGN are picked up using this technique but they are much less common. 

While our sample is small and we are deahng with small number statistics, our results 
confirm several trends observed before: 1) ULIGs are almost invariably mergers or interacting 
galaxies; 2) the fraction of LIGs undergoing mergers is significantly less; 3) objects with 
higher FIR luminosity are more likely to contain AGN. 

We find, on the other hand, a larger fraction of exponential or near-exponential profiles 
(nearly half of the sample) than in previous surveys. As discussed in Section 4.1, many of 
these surveys have based their classification of profiles on one dimensional fits which can be 
subject to errors due to isophote twists or other small perturbations. It is possible that the 
true fraction of objects with exponential profiles may be larger, as suggested by our study. 
However, because of the small size of our sample this is simply speculation at this stage. 

We do confirm, however, that modeling in two dimensions is very effective in highlighting 
features that could be easily missed by visual inspection or one dimensional fitting. Residual 
images reveal details such as double nuclei, dust lanes, tidal debris, and secondary cores that 
allow us to identify more accurately those objects that are perturbed. 

Our results show the effectiveness of using adaptive optics systems in combination with 
two dimensional modeling to study morphologies of infrared galaxies. Future morphological 
surveys of this kind can be done with success for other types of interesting galaxies. 
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Keck Foundation. This work was supported in part under the auspices of the U.S. Depart- 
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Fig. 2. — Lick and Keck near IR AO images of the central regions of each galaxy in the 
sample. 
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Fig. 3. — A BPT diagram of emission line flux ratio [0III]/H/9 versus the ratio [NII]/Hq;. 
From left to right the galaxies plotted are: (1) FF 1656+2644, (2) FF 1110+3130, 
(3) FF 1709+5220, (4) FF 0934+4706, (5) FF 1712+3205, (6) FF 1721+2951, 
(7) FF 1517+2800, (8) FF 1138+4405, (9) FF 1723+3845, (10) FF 1318+3250, 
(11) FF 1412+4355, (12) FF 1429+3146, (13) FF 0834+4831, (14) FF 1519+3520, 
(15) FF 1651+3001, (16) FF 1122+4315, (17) FF 0835+3142. 
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Fig. 4.— Optical spectrum of the z = 0.2613 ULIG FF 0819+2707 in rest frame. 




Fig. 5.— Panel a is an image of the z = 0.2630 ULIG FF 1708+4630 before subtraction. 
Panel b is an image of the object after subtraction, showing a secondary nucleus or merging 
companion. 
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Fig. 6.— Residual image (model subtracted) of the z = 0.1484 LIG FF 1122+4315 showing 
much structure near the AGN nucleus. 
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Fig. 7. — Two merging disk galaxies form the z = 0.1353 LIG FF 1412+4355. The image 
has been smoothed using a gaussian with a — \ pixel to highlight the tidal tail and debris. 
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Fig. 8.— Keck NIRC-2 Narrow Camera image of the z = 0.1806 LIG FF 1429+3146, dis- 
played in a log scale. Panel a shows the galaxy before subtraction and panel b after sub- 
traction. The dark regions are an artifact of over-subtraction due to the highly perturbed 
morphology and multiple components of this galaxy. 



